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Abstract 

We discuss gluon production by the Schwinger mechanism in colhnear color- 
electric and magnetic fields which may be realized in pre-equilibrium stages of ultra- 
relativistic heavy-ion collisions. Fluctuations of non-Abelian gauge fields around a 

■ purely color-magnetic field contain exponentially growing unstable modes in a lon- 
D . gitudinally soft momentum region, which is known as the Nielsen-Olesen instability. 

i ' With a color-electric field imposed parallelly to the color-magnetic field, we can for- 

mulate this instability as the Schwinger mechanism. This is because soft unstable 
, modes are accelerated by the electric fields to escape from the instability condition. 

I Effects of instability remain in the transverse spectrum of particle modes, leading 

■ to an anomalously intense Schwinger particle production. 

1 Introduction 

Multi-particle production in strong fields is a typical unstable phenomenon that can be 
seen in extreme situations such as ultra-relativistic heavy-ion collisions, high-energy as- 
^ ■ trophysical objects, and possibly high-intensity laser. This phenomenon is in general 

nonperturbative even for weakly interacting systems because the coupling with external 
strong fields compensates the weakness of the interaction. A straightforward description 
of such an instability may be given by the effective action of the strong fields. A famous 
example is the Euler-Heisenberg action, which is an effective action of strong electromag- 
netic fields due to electron's one-loop contribution. However, such an effective action 
does not directly tell us the information of dynamical processes which cause instability. 
To obtain dynamical pictures of instabilities, we have to resort to other methods which 
allow us to describe the degrees of freedom relevant for instabilities. 

When the produced particles are normalizable and asymptotically stable modes, one 
can describe particle production within the canonical quantization. This is indeed the case 
with the Schwinger mechanism [1] where a pair of a charged particle and an antiparticle 
is created in a strong field. 
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On the other hand, when a system contains genuine unstable modes which are unnor- 
malizable and keep growing even at asymptotic regions, one cannot apply the canonical 
quantization to those modes. When the mode amplitude grows in an unnormalizable way, 
this mode should be treated as not a particle but a field. However, if the instability lasts 
only for a finite duration, we have, after the instability, an asymptotic region where stable 
particle modes can be defined. Then we can provide particle interpretation of the modes 
by the canonical quantization, and therefore the instability can be viewed as a particle 
production phenomena. In fact, this situation is seen in the parametric resonance, where 
the time- dependent external field induces instability only for a finite time interval, and 
particle number is defined in the free regions [21 13]. 

Notice that one encounters both of the cases in non-Abelian Yang-Mills theories; The 
one-loop effective action of non-Abelian gauge fields under constant electric or magnetic 
background has an imaginary part, which indicates the existence of instabilities [1]. The 
instability in a purely electric field consists of normalizable and asymptotically stable 
modes, and thereby can be interpreted as particle productions O |6]. One can formulate 
it as a non-Abelian analog of the Schwinger mechanism [7]. In contrast, the instability 
in a purely magnetic background is caused by rapid growth of particular fiuctuations and 
thus cannot be formulated as particle production unless the magnetic field is turned off 
in finite time so that the modes are stabilized. This "Nielsen-Olesen (N-0) instability" 
[8] is characteristic of non-Abelian gauge fields because the properties of self-interaction 
and spin-1 are necessary. 

The instabilities in these two limiting cases have been individually studied in the 
context of early time dynamics of ultra-relativistic heavy-ion collisions [H (TUl IH [12] (see 
also Ref. [13] )• However, the strong non-Abelain fields that appear in reality (called a 
'glasma') would consist of longitudinally extended fiux tubes in which both the electric and 
magnetic fields are non- vanishing and collinear along the beam axis p3]. It is important 
to investigate particle production in this configuration to understand time evolution of 
created matter towards a thermalized quark-gluon plasma. In this Letter, we study, within 
the canonical quantization approach, the gauge field instability (gluon production) when 
both electric and magnetic background fields are present. 



2 Formalism 

We consider the SU(A^) pure Yang-Mills theory, and decompose the gauge field (a = 
1, ■ ■ ■ , N"^ — 1) into a classical background and a quantum fiuctuation around 
the background: = -|- We take a covariantly-constant background field [T5] 
satisfying 

D^F'^^'' = Q, (1) 

where is a covariant derivative with respect to the background: D^x"" = {df^S"''^ — gf"''"^^''^) y^-, 
and is a classical field strength: = dy^A^ - d^A^^ - gf^^A^^Al. We employ the 
background covariant gauge [7]: 

C = -\f;,F''^- - ^ [D.A'^Y - I {D'^C'') D^C% (2) 
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where C is a gauge parameter (to be set as C = 1) and C"" (C"^) is the Faddeev-Popov 
(anti-)ghost. The last covariant derivative Z?^ contains both the background and the 
fluctuation In order to be consistent with the one- loop calculation of the effective 
action, we retain in the Lagrangian up to quadratic terms with respect to the quantum 
fields. Then, the equations of motion for the quantum fields are linearized as follows: 

D^D'^A''^ - 2gr^''F^>"'Al = , (3) 
D^D^'C = , D^D^'C" = . (4) 

We consider the classical background field which gives spatially constant field strength. 
Then, one can express the field strength as [151 E] 

F;, = F,,n\ (5) 

where F^i, is an Abelian field strength and is a color vector, normalized as n°'n°' = 1. 
The field strength ([5]) is given by the gauge potential 

K = A,n\ (6) 

where is an Abelian gauge field giving F^j, = dfj^Ay — d^A^. This is actuafly a solution 
to Eq. ([T]). The global residual gauge symmetry allows us to rotate the color vector n°- 
into the Cartan subspace of SU(A^): Un°'T"'W = h-^H^, where f/ is a constant SU(A^) 
matrix, are generators of SU(A^) and (A = 1,2, ■ ■ ■ , N — 1) are elements of the 
Cartan subspace: G {T"| [T",T^] = 0}. 

We expand the SU(iV) space by H"^ and v^T"", instead of T". Here, are eigenvectors 
of ad{if^}^'^ (i.e., H"^ in the adjoint representation): 

ad{i7^}''X = • (7) 

The corresponding eigenvalues a"^ may be regarded as an (A^ — l)-dimensional vector, 
which is called the root vector. According to this Cartan decomposition, we re-organize 
the quantum fields A'^, and (7° as 

{A-,} ^ {a^} ® {W^ ^ v^A-^} , (8) 
{C"^}^{c^}©{ry" = ^-C"}, (9) 
{C"^} ^ {c^} © {ry° = v^^C^} . (10) 

This decomposition simplifies the interaction between the fiuctuations A^^, C"", (7" and 
the classical field ^4°. The fields in the Cartan subspace, (p^ = a^, c^, c^, are free from 
the background field: D^j^cj)^ = d^cj)^, and thus we need not consider these fields in the 
following. For the off-diagonal fields, = Wj^ ,7]°" ,7]°" , the covariant derivative reduces 
to an Abelian form as 

5^$"= (a^ + ze^A^) <l>^ (11) 

with Cq, = —n^a^^g. Gauge invariance of has been exphcitly verified for SU(3) [16]. 
One can discuss quark production in a similar way [T7]. 

In terms of these new fields, the equations of motion ([3]), (jl]) are rewritten as 

{d, + lec^A^fw^ + 2iec,F^,W^ = , (12) 
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{d^ + ieaA^) r/" = 0, {d^ + ieaA^) ^ = 0. (13) 
Solutions of these equations may be expanded by normal modes: 

cr=±,L,S n 

r{x) = j2[hn{x)c: + h:{x)d"j], (m) 

n 

v'^ix) = J2[hnixK + K{x)d:^], (15) 

n 

where n is quantum number(s) characterizing the mode, such as momentum p, and a{= 
±,L,S) denotes polarization. fn\x)e'l^'^ and hn{x) are c-number solutions of the equations 
of motion (IT^ and (1131) . respectively. e^T'' are polarization vectors satisfying the orthonor- 
mal condition e|r^*e^'^-'^ = —fj'^'^ and the completeness condition X^o-t ^a*^^^'^'^^^'^^* ~ ~9iJ.v 
with the metric ff'^'^ such that f]^^ = fi = ff"^ = f]^^ = 1 and other components are 
vanishing. We identify the solutions fn{x)e^i^^ and hn{x) as positive frequency modes, 
and f!f^*{x)e^^^ and h^{x) as negative frequency modes. The positive frequency solutions 
are to satisfy the following normalization conditions: 

where the inner product is defined as 

(01, h) = ^ j d^x {(PI ■ Doh - Do^l ■ 02) . (18) 

Notice that the norms of the negative-frequency modes have the sign opposite to those 
of the positive-frequency modes, and the positive and negative frequency modes are or- 
thogonal to each other. With these orthonormal conditions, the canonical commutation 
relations imposed on new fields VF^,r7°,f/" yield the following commutation relations: 

[a(^)^ a(;)^t] = 5M/3t] = (19) 

{C,c^J} = {d':^,d^J}=t6^p6nm. (20) 

Now we can regard the operators a^T''", fe^T''", c", d", and J" as annihilation operators 
for particles or antiparticles. 

3 Instability and gluon production 

With the canonical framework defined above, we investigate the quantum dynamics of 
fluctuations in the coUinear color-electric and magnetic background field: 

E" = (0, 0, E)n^ B" = (0, 0, 5)n", (21) 
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which is reahzed by A^^ = {0,—By,0,—Et). It is convenient to choose the polarization 
vectors as eigenvectors of the field strength F'^'^. Then, all the equations that the mode 
functions fn (x) and /i„(x) follow can be put together into a Klein-Gordon- like equation 




{d^ + ieaA^Y + <l>{t,x;m'^) = 0, (22) 

where $ = In the following, the subscript a of Cq, will be often omitted. The 

"mass squared" represents the effects of spin-electromagnetic field interaction and is 
given by 

for gluon modes with a = ± 
2 _ I 2ieE for gluon modes with a = L , . 

' o„-^ z? ^Qj- gluon modes with a = S 
for ghost modes. 

Because the a = L,S modes acquire pure imaginary m^, we should suppose that the 
negative frequency solution for a = L mode is fn , and vice versa. 

Let us first consider the case where we have only a color-electric field {E ^ 0, B = 0). 
Ambj0rn and Hughes [7] have studied the Schwinger particle production in the canonical 
quantization approach in the background covariant gauge ([2]) . One can find exact solutions 
of Eq. (l22l) and show that all the mode functions satisfy the normalization conditions 
(fT6|) and (ITTl) . Therefore, in the pure electric field the canonical quantization can be 
completed to give particle interpretation of the fields. However, under electric fields, 
selection of the positive and negative frequencies is not unique. There are infinite numbers 
of solutions satisfying the conditions f|T6l) and f|T71) . If there is no electric field, we can 
choose the positive and negative frequencies uniquely as e^*'^''* respecting the translational 
invariance in time. In contrast, under electric fields, the system loses that invariance, so 
that we cannot select a specific set of solutions without other criterion. To decide proper 
solutions, asymptotic WKB criterion is employed [3 [TIE]. In this criterion, the positive 
and negative frequency solutions are selected so that they behave at asymptotic region, 
|t| — )• oo, as ^(t,x) ~ exp(zS'cz) with Sd being the classical action of a charged particle 
under the electric field. As a direct consequence of particle production, positive and 
negative frequency solutions defined at t — > — oo and those at t — > +oo are different. The 
former, which are referred to as in-solutions, are 

-fax „ipx 

with a± = m]_ = rn? +pl+Py and ^ = w ^(p^ + cEt). D^{x) is a parabolic cylinder 
function. The latter, out-solutions, are 

-fax pipx 

^ (2eE)3 ^ ' ^ ^ 7(2^ 

Both solutions have the structure of plane waves distorted by the electric field. Associ- 
ated with these two kinds of mode functions, we get two kinds of particle creation and 
annihilation operators, ajT'*™ and alp^°^^, etc. Accordingly, two kinds of vacua |0,in) and 
|0,out) are defined. The creation and annihilation operators of in and out are related 
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with each other via Bogohubov transformations. Then one can relate the in-vacuum 
and the out- vacuum and calculate the vacuum persistence probability |(0,out|0,in)p. 
It has been shown that contributions from the unphysical modes, a = L,S, and those 
from the ghost modes are totally canceled out, and the pair creation rate w defined by 
1(0, out|0, in) p = exp{—VTw) with V and T being space and time volume, respectively, 
is just twice as large as that of a massless charged scalar field [7]. Also, we can calculate 
the momentum distribution of out-particles condensed in the in-vacuum, in other words, 
the momentum distribution of particles which are created during infinite-time imposition 
of the electric field. For the physical modes, a = ±, we find 

(0,in|aM"-tt^H-ut|Q^i^^^^^p^_^^_^ (26) 

which depends only on the transverse momentum p]_ = p^. +^^0 This result is the same 
as that of a massless charged scalar field. In contrast, for the unphysical modes, a = L,S, 
and the ghost modes, the expectation is vanishing because of their unusual commutation 
relations Ipl) and fl20l). 

In contradistinction to the electric field, the magnetic field has no dynamical effect. It 
just discretizes the transverse momenta into the Landau levels as 

pl^{2n + l)eB (n = 0, 1, 2, ■ ■ ■ )• (27) 
Under the pure magnetic field {E = 0, B 0), one can find a solution of Eq. ( 122|) as 

%{t, X- w") = ^=e---Vn(2/)7^e^^^^+^P^^ (28) 
■s/2ujp 27r 

where ujp = ^ym"^ + pl + {2n + l)eB and v'„(?/) oc Dn (^^^^{eBy + p^)^ . The spin- 
magnetic field interaction ( !23l) (with E = 0) shifts the Landau levels as shown in Fig. [H 
Notice that the transverse mass squared defined as m]_ = m? + {2n + l)eB is negative 
for the lowest Landau level n = with a = + mode, and the one-particle energy reads 
ujp^^n=o = \/pl — eB. This is pure imaginary when \pz\ < y/eB. Therefore, the mode 

functions show exponential growth ~ exp (^y/eB — p^t^- This is the N-0 instability |8]. 

When this instability occurs, the normalization ([T6|) cannot hold, so that we can no longer 
gain particle interpretation of the field. If the magnetic field is not constant in time and if 
there are asymptotically free regions where particle modes can be defined by the canonical 
quantization, the N-0 instability could be described as particle production phenomena. 
However, as far as the authors know, there has been no explicit demonstration of it. We 
will show that we can define asymptotically stable particle modes even in the presence 
of the constant magnetic field provided there is a constant electric field parallel to the 
magnetic field. 

Notice that the N-0 instability occurs only at low longitudinal momentum region 
\pz\ < V eB. Under the collinear electric and magnetic fields {E 0,B 0), particles are 
accelerated by the electric field and their momenta increase in time as Pz = eEt + const. 



*^As discussed in Ref. [5], Pz-dependence of the distribution cannot be obtained by the method of the 
asymptotic WKB criterion. 
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Figure 1: A schematic of the Landau level and the level splitting due to spin-magnetic 
field interaction 



Therefore, even if a particle lies in the unstable momentum region \pz\ < \/ eB at some 
time, it escapes from there unfailingly at later time. Owing to this mechanism, the expo- 
nential growth of the mode function is regularized in the presence of the color-electric field, 
and we can carry out the canonical quantization and obtain the particle interpretation 
for all the modes. Indeed, exact solutions of Eq. ( |22l) under the coUinear fields 



^'^{t, X- n?) = -—^D,,^_r (-e-i-e)yn(l/) , (29) 

$r ^; m') = j^D_,a^-^^ {en)My)^^, (30) 

with a„ = + I) ^ do not diverge at t — j- +oo and satisfy the normalization conditions 
(ITB]) and (IT7I) . This is true even for the lowest Landau mode with n = and a = +. 
Since the magnetic field has no dynamic effect, the selection of the in- and out- mode 
functions is not affected by it. Therefore, the Bogoliubov relations are almost the same 
as those in the pure electric field. The differences are the emergence of the Landau levels 
fl27|) and the level splitting between the two transverse modes (see Eq. fl23l) ). Because 
for the unphysical modes are independent of B, cancellation between those modes 
again holds even in the presence of the magnetic field. Actually, the pair creation rate 
has contributions only from the two physical transverse modes: 



5^5^1n[l + e--(2n+i-2<x)fj _ (3^) 



2 (2tt) 

^ ' CT=± n=0 

Furthermore, the in-vacuum expectations of out-number operators for the unphysical 
modes are all vanishing also in the presence of the magnetic field. Those for the physical 
transverse modes are 

(0, in|aW--*t«(±)«out|Q^ .^^ ^ ^"^^'"""'^'^^ (^- ^32) 

In Eqs. (pTl) and (1321) . the lowest level (a = + and = 0) is worthy of great remark. 
In this mode, the effective mass square is negative, mj_ = —eB, so that the index of the 
exponential factor is positive: 
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This means that an anomalously large number of gluons are produced in the lowest 
mode when compared with the usual Schwinger mechanism of QED or scalar QED. This 
is because the field fluctuations amplified by the N-0 instability are converted to real 
particles by the electric field. 

The coUinear color electromagnetic field fl2T|) may be realized in flux tubes created at 
the initial stage of heayy-ion collisions The strongly enhanced particle production 
causes rapid decay of those coherent fields into particle degrees of freedom, which would 
be followed by the formation of a quark-gluon plasma. 



4 Discussion 

Our analysis has relied on the linear approximation. If the field amplification by the N-0 
instability continues for a long time, the linear approximation would get invalid. A typical 
time scale of the N-0 instability is given by the inverse of the growth rate ^no ~ l/V eB. 
Therefore, the time scale when the non-linear corrections become important would be 
^ni ^ 1/ VeB. Actually, it has been demonstrated by numerical calculations that the 
exponential growth of the gauge fluctuations stops at some time (later than tNo) due to 
non-linear corrections |18]. Meanwhile, the time scale of the stabilization process by the 
electric field teie is given from the condition that a particle is accelerated away from the 
instability condition Ip^l < \/eB, so that teie ~ VeB/eE. If the electric and magnetic 
fields are the same order of magnitude, which is a situation expected in a glasma, teie is 
the same order as t-^o, so that the electric field plays a dominant role to stabilize the N-0 
instability. 

To understand more clearly our result that the N-0 instability is stabilized by the 
electric field and induces the strong particle production, the following simple toy model 
which describes particle production by an instability may be useful. Let us suppose that, 
without specification of background fields, a real scalar field suffers an instability during 
a finite period, < t < T. The mode functions of the field have the following properties: 

(i) for t < the modes are stable and have the plane wave solutions; (i^ > 0), 

(ii) for < t < T the modes are unstable and show exponential growth or decrease; e^'*'* 
(7 > 0), and (iii) for T < t the modes are again stable; ^e=F^^* (fi > 0). The canonical 
quantization can be performed at t < and at t > T, and correspondingly two kinds 
of particle definitions of 'in' and 'out' are obtained. One can find the relation between 
these two kinds of particle definitions by smoothly connecting the mode solutions in the 
three time regimes. If the solution -^e"*'^* for t < is smoothly connected to the lin- 
ear combination of the solutions for < t < T, and is subsequently connected to those 
for t > T as -^^^ (ae"**^* + /3*e*™) , then the particle annihilation (creation) operators of 

in-particle, related with those of out-particle, a°^^^'^\ by the Bogoliubov transfor- 

mation: a°"* = aa'° -|- /9a'°^. One can explicitly confirm the condition |ap — = 1, 
which guarantees the unitarity, holds irrespective of the values of u, Q, 7 and T. The 
number of particles created by the instability is 

(O,in|a-*t«out|o ^ = ^ ^(^ + ^](^ + 7^] ^"'^ ■ (34) 

T>l/7 16uj \ 7 / V " / 

Thus, the number of created particles increases exponentially with respect to the time 
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interval of the unstable region, T. In this simple demonstration, the existence of the stable 
modes in asymptotic regions has been crucial to formulate the particle production in a 
canonical way. If we compare this simple model with our original problem of the gluon 
production, the instability parameter 7 corresponds to 1/tNO ~ \J eB. Though there is 
no explicitly stable region in temporary direction under the constant magnetic field, the 
electric field provides us with the stable regions in the momentum space, i.e. \pz\ > VeB- 
Hence, T should be replaced by teie ~ VeB/eE. By substituting these parameters, we 
find the exponential factor e^^'^ in Eq. roughly reproduces the exact result ( l33l) . 
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